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Abstract. We provide a rigorous microscopic derivation of the thermodynamic descrip- 
tion of equilibrium crystal shapes in the presence of a substrate, first studied by Win- 
terbottom. We consider finite range ferromagnetic fsing models with pair interactions 
in dimensions greater or equal to 3, and model the substrate by a finite-range boundary 
magnetic field acting on the spins close to the bottom wall of the box. 
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1. Introduction and Results 



1.1. Winterbottom variational problem. The Winterbottom theory [ Wi] gives a phe- 
nomenological prediction for the equilibrated shapes of small (that is disregarding grav- 
itation) crystal particles placed on solid substrates. The notion of equilibrium comes on 
two levels: 

(i) As far as the bulk properties are considered, both the particle and the vapor around 
it are assumed to represent two different thermodynamic phases of the same physical 
substance. 

(ii) The total interfacial free energy of the system comprises the anisotropic surface tension 
between the particle and the vapor, as well as the boundary surface tensions between the 
particle and the substrate and, respectively, between the vapor and the substrate. The 
substrate could, thus, exhibit a preference towards one of the phases, and the macroscopic 
Winterbottom equilibrium shape corresponds to the state of minimal interfacial energy at 
a given particle volume. 

We use (5 to denote the inverse temperature at which the two phases (crystal and vapor) 
coexist, and the parameter rj to quantify the preference of the substrate towards one of 
the two pure phases. In the magnetic interpretation of the microscopic models we shall 
consider below r\ is the boundary magnetic field. According to the suppositions (i),(h) 
above, the two principal players to determine the equilibrium shape are the anisotropic 
crystal- vapor surface tension rp = r^cv : i— * and the number = r^sc - r /3,SV) 
where r^sc and r^sv are, respectively, the substrate-crystal and the substrate-vapor 
boundary surface tensions. Indeed, if the crystal particle occupies the region P of M. d 
its boundary splits into two disjoint pieces dP = dPsc U 9Pcv (Figure ||), and, up to a 
constant, the total surface energy of the system can be written (at least for piece-wise 
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Figure 1. Left: The crystal P on the substrate. The boundary is decomposed into two 
pieces, corresponding to the crystal/vapor interface and the crystal/substrate interface. 
Right: The Winterbottom shape Kp tn (shaded region) is obtained by intersecting the 
Wulff shape and the half-plane; on the picture A/3 >r; is negative. 

smooth shapes) as 

W^{P) = [ Tp(n x )dS(x) + Ap jV \dPsc\ ■ (1-1) 

JdPcv 

The equilibrium particle shapes at a prescribed volume v correspond, in the above nota- 
tion, to the minimizers of the following variational problem: 

(WP)^ Wp >T , (P) — > min Given vobj (P) = v. 



Winterbottom [Wi| described the geometric construction of the solutions to (WP)^ 



which, in fact, is an easy consequence of the general Wulff construction [Wu| of the equi- 
librium crystal shapes: First of all the minimizers are scale invariant under dilatation. In 
order to construct the unnormalized Winterbottom shape pick the "free" Wulff crystal 



Kp = p| {x£R d I (x,n) <T,j(n)} 



and intersect it with the half-space 

JCp, v = /C^pljxer* I x d > - A^} , (1.2) 

where we use the coordinate representation x = (x\, Xd) of a vector x G M. d . It is 
convenient to parameterize the substrate surface by one of the coordinate hyperplanes, 

say by V 
given by 

■ V \l/d 



say by Va = {x : Xd = 0}. The equilibrium Winterbottom shape at volume v is, then, 



A / v\ L / a 

KpM = [j) (Afrtied + lCp,n) , (1-3) 

where is the unit coordinate vector in the d-direction and A = \((3, rj) = void (]Cp >r) ) is 
used to denote the volume of the unnormalized Winterbottom shape. 

As we shall explain in Subsection |l.5| , the value of the difference A^ jr) always lies in the 
interval [— Tp(ed), 7/3 (e^)]]]. Thus, one of the following three cases happens: 

1. Complete drying: Ap^ = rp(ed)- In this case Kp tTj = ICp, which means that the 
Winterbottom droplet has no energetic reasons to be near the substrate surface and 
can be located anywhere above it. 

2. Partial drying /wetting: In this regime Ap >r) 6 (— T^(e^), T^(e^)). Accordingly, the 
convex body KLp^ is obtained by removing a proper cap from Kp (Figure |lj). We 



1 We shall always assume that Tp(x) = t@(— x). 
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shall, furthermore, distinguish between the partial drying, Ag „ > 0, and the partial 
wetting, Af3 jTI ^ 0, cases. 
3. Complete wetting: = — Tg(e*d). In this case, /Cg jr) has volume zero and, therefore, 
cannot be rescaled to get a volume v. Physically this corresponds to the formation 
of a microscopic film along the substrate, preventing the vapor phase to reach it. 

Our main objective here is to give a rigorous derivation of the phenomenological Winter- 
bottom picture in the scaling limit of microscopic models of lattice gases, that is directly 
from the basic principles of the Equilibrium Statistical Mechanics. The class of microscopic 
models we consider are ferromagnetic Ising type models with finite range of interactions. 
In the subsequent subsections we shall briefly recall these models as well as the ther- 
modynamic procedures leading to the notion of bulk phases and various surface tensions 
which appear in the functional Ql.l| ). The results fall in the framework of the Li-theory 
of phase segregation, which has been initiated in the worksJPl on the nearest neighbor 



Ising model and in the works jABCPfl , [|BCP|| , flBBBPl , flBBPfl on the Kac Ising models. 
Recent developments of the theory have been prompted by the work Ce] and comprise 
a rigorous derivation of the (free) multi-dimensional Wulff crystal shapes in the context 
of Bernoulli bond percolation |Ce{| and in the context of nearest neighbor Ising model in 



[Bl] and | CePil | . Most recently the theory has been generalized to the case of symmetric 



^-states Potts models [|CePi2| . We refer to the review article fBIVlf for a detailed account 
of the Li-theory and its relation to the Dobrushin-Kotecky-Shlosman approach (see e.g. 
[ DKSl , Q, 0, (I|, gVg , gScjl) to the problems of phase segregation. 



1.2. Other microscopic derivations of the Winterbottom construction. The deriva- 
tion of the macroscopic geometry of an equilibrium droplet in contact with a substrate 
has been accomplished in several models, which can be divided into two classes. 

The first class is that of SOS-type models, with a positivity constraint, a fixed volume 
under field, and an attraction to the wall. The first such study was done for the 1 + 1- 
dimensional model in JCDRJ (see also |MR|] for an alternative approach), while its higher 
dimensional counterpart was treated injjBj]; see also [ CMR| for an analysis of the effect 
of substrate roughness in 1+1 dimensions. Other related works, which neglect the inter- 
action between the substrate and the interface (therefore describing a completely neutral 
substrate) , are flBD|| , DG1 ] and [ DM ] ; notice however that these works cover only one spe- 
cial point of the full regime of complete wetting. This class of models provide a somehow 
appealing description of this phenomenon, but has several shortcomings. The first one 
is that the full phenomenology of the wetting transition cannot be observed, due to the 
very particular constraint that the droplet has to be the graph of a function (e.g. it is of 
course impossible for the droplet to detach from the wall). The corresponding variational 
problem also displays some non-physical properties, e.g. the minimizer are in general not 
scale invariant. And of course, these SOS models are effective models for interfaces, not 
truly microscopic ones, so they do not clarify the mechanism which generates the inter- 
face itself. On the other hand, due to these simplifying features, it is possible to obtain 



some additional information about fluctuations around the equilibrium shape, see [DGI], 
which go beyond what we can achieve for the models considered in the present work. For 
an in-depth review and references on these effective models, we refer the reader to the 
survey |G| . 

The second class contains more realistic lattice gas models of the type we consider here. 
Only the two-dimensional nearest-neighbor Ising model had been studied before [PV1 



PV2| ], but the description provided is rather complete. The wetting transition can be fully 
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described, and the characterization of typical configurations is very accurate, especially 



when these results are combined with the local limit estimates of [ISc]: There is a unique 
macroscopic droplet whose boundary is close in Hausdorff distance to one of the solutions 
of the Winterbottom variational problem, and all the other droplets are of size at most 
log L, L being the linear size of the box. Notice that at low temperature, it is also possible 



to obtain precise informations on the fluctuations around the minimizer, see [DH 



1.3. Microscopic model: Bulk phases. We start by introducing the measures describ- 
ing finite-range ferromagnetic Ising models with pair interactions in the bulk; the complete 
description, including the influence of the substrate, is given in Subsection |l.5| . 
The system is contained in the finite box 

Bjy = {* G Z d : \i k \ ^N,k = l,...,d}. 

The interactions are given by a set of real numbers J = (Jjj)i j£z d > with the following 
properties: 

^ T ^ T 

2. the graph (Z d , Ef) with edge set £j consisting of all (unoriented) pairs of vertices 

with J\i-j\ > is connected; 

3. J fc ^ for all k G Z d ; 

4. J/% = if \\k\\i > R, for some R < oo (the smallest such R is called the range of the 
interaction) . 

Let a G {—1, l}^ ■ The Hamiltonian in Bjv with boundary condition a is given by 
H CT (cr) = -- ^2 JijOiOj- ^ Jijvi&j, 

{i,j)CB N i6Bjv,i?"' 



and the corresponding (finite-volume) Gibbs measure at inverse temperature (5 is the 
probability measure on {-1,1} Z defined by 



exp (-/3H ct ((t)) if <Xj = <7i, for alH 
I otherwise. 



6 N, 



Two boundary conditions of particular importance are the + b.c. defined by W = 1, 
and the — b.c. a = —1; the corresponding finite-volume Gibbs states are denoted by 
H N , and ^P N _ respectively. It follows from correlation inequalities that the weak limits 

/x+ = limjv^oo^^ + and = \\m^^ 00 ^ N _ exist and are translation invariant (in fact 
ergodic). It is a basic result of Statistical Mechanics that in dimensions d ^ 2 there exists 
< (3 c (d, J) < oo such that yP+ = fi_ if (3 < (3 C (this also implies uniqueness of the infinite- 
volume Gibbs state), while [i^ ^ if > j3 c - The order parameter associated to this 
phase transition is the spontaneous magnetization, 

m*(/3) = lim //J I — !— V <n J = /xJ(<t ) • 
\ 1 1 ieB N J 

It can indeed be shown that m*((3) = when [3 < f3 c and m*{(3) > when (3 > f3 c . Of 
course, by symmetry, ^_{oq) = —m*{j3). 
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1.4. Microscopic model: Surface tension. Let ft 6 § d 1 be a unit normal. Assume 
for the definiteness that (n, e^) > \f\f~d. Given two positive real numbers L and M, define 



a j |(x,e fc )KL/2 for k = l,...,d-l\ 

K M , L (n) = (^«" and|(n,,)KM/2 ^ 



Thus, AM,i(n) is a parallelepiped with the base orthogonal to n and having area L d 1 /(n, e^) 
and height M. With a slight abuse of notation we identify A^i(n) with its intersection 
with Z d ; Am,l(^) = ^-M,L(n) nZ d . Let ^X. M L (fC\ anc ^ ^A lfi (n) be the partition functions 
on Am,l(tj) with respectively "+" and n boundary conditions, the latter being defined by 
(Tj = sign((n, i)), with sign(O) = 1. The (per unit area) surface tension of the it-interface 
stretched in the direction orthogonal to n is defined via: 

r,(n) = - L hm log ^— (1.5) 

The important ferromagnetic feature of the model (or, possibly, rather of the proof of the 



corresponding fact) is that the limit in (1.5) is defined and does not depend on the order in 



which the numbers L and M go to infinity. In [MMR] this result has been formulated for 
the sequences of domains tending to oo in the sense of Landford, but the corresponding 
proof goes through also in the limM^oo lim^oo case. 

Finally, as it has been proven in [LePi], the surface tension is uniformly strictly positive 
whenever (3 > f3 c . 

1.5. Microscopic model: Wall free energy. Microscopic models for the substrate 



surface have been studied in [ FP1 ] and [ FP2 | in the context of the nearest neighbor model. 



Let us reformulate, with brief comments on the validity in the general case we consider 
here, all those of the results of [ FP2| which we are going to use in the sequel. 



Define the lattice half-space h d = {i S 7L d : > 0}. Our spin configurations a are the 
elements of {-1,1} L . The microscopic influence of the substrate is modeled by magnetic 
fields (chemical potentials in the lattice gas language) 77 = (771, ... ,r/ r ) acting on spins in 
the first r microscopic layersf^of L rf . Thus, the formal Hamiltonian on {— 1, 1} L is given 
by: 

r 

IF((7) - - ^ JijOiUj - ^2 ^2 rj k ai. (1.6) 

(i,j)Ch d k=l i£h d :i d =k 

Given two different vectors r\ and 7/ of boundary magnetic fields we say that 77 is larger 
than 7/ if 77^ ^ rj' k for all k = 1, . . . , r. 

In order to describe the appropriate surface energy corrections induced by the substrate, 
we put it in contact with either of the two "— " or "+" bulk phases. This is done in a 
standard way: Given A £ N we construct the finite lattice box 



{i £ L rf : -N < i k < AT, k = 1, . . . , d - 1, < i d ^ N} , 



and consider the models with formal Hamiltonian (|1.6| ) and, respectively, with "— " and 
"+" boundary conditions on L d \ B/y. Let //^ _ and /jp^ + be the corresponding finite 



2 We do not assume any connection between the number r and the range of the interaction {Jij}. 
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volume Gibbs states on {—1, 1}° N and _, + be the associated partition functions. 
The difference A^n between the interfacial free energies is defined via: 

1 Z P ' V 

A«„ = lim . —7-3 : log — a N ' + . (1.7) 
Pfl Af-^oo (2iV) d - 1 



Proposition 1.1 (| FP2j ). The limit in (1.7) is we// defined and monotone non-increasing 
in r]. Furthermore, uniformly in r] 



I A/3,,, | < T^(e rf ). 



li 



As in [ FP2 | the assertion of the theorem follows from the FKG properties of the ferro- 
magnetic measures. We refer to Section 3.1 in FP2fl (with obvious modifications to fit the 
general case we consider here) for details. 

Remark 1.1. On 

nario under the (J^ 
along the substrate 
definition of Tp{ed) 
values of rj, which 
Gibbs state |FP2fl . 
rj) is, however, an 
results for effective 



the heuristic level the bound in Q1.8D should be clear: one possible sce- 
^ _ measure is to create a microscopically divergent film of "+ " phase 
surface reducing, thus, the ratio in (|l.7| ) to the one appearing in the 
in (|l.5|) . This scenario certainly becomes dominant at large positive 
happens to be equivalent to the uniqueness of the limiting boundary 
A non-trivial analysis of the corresponding phase diagram (in terms of 
almost entirely open question, c.f. |Ch|, | BD2f , |CVf| for the related 
interface models. 



1.6. The result. In principle, the results of the paper should hold for any (3 > (3 C . 
However, our approach relies on the validity of Pisztora's renormalization machinery, which 
we shall describe in detail in Section ||. Consequently, we need the following assumption 
on the inverse temperature (3: 

(A) Pisztora's coarse graining holds at the inverse temperature (3. 

We use 25 to denote the set of (3 > (3 C for which the above assumption holds. 

Define the total average magnetization on D^r as 



Mi 



i€D N 



a;. 



The averaged magnetization Mu^ concentrates (as N — > oo) under the measure //p + 
on the spontaneous "+" phase magnetization m*{(3). In other words the bulk phase 
induced by ^ v + is the purely "+" one. In order to enforce phase segregation one should 
introduce a canonical type constraint which would shift the average magnetization Md n 
inside the phase coexistence region. In the lattice gas language this would amount to fixing 
the average number of particles strictly between the two extremal equilibrium densities. 
The canonical constraints we are working with here are of the integral type: Given m £ 
(— m*(f3), m*{(3)), let us define the canonical conditional measure 



^K,m,+ ( CJ ) 



1 



J Bjv,m,- 



■ cxp 



1 {M ]0n ^m}(°) 
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The impact of such a conditioning on the bulk properties should amount to the creation 
of a "— " phase island of macroscopic size close to 

. . a m* — m , 

v (m) = — . 1.10 

v ; 2m* y ' 

Our microscopic justification of the Winterbottom construction above the complete wet- 
ting threshold; A^ > — Tp(ed), gives a rigorous meaning to the heuristic picture above 
and, moreover, asserts that under the natural scaling the shape of the "— " phase island 
converges to the Winterbottom shape JC/3 >V (v(m)). Such a result is, of course, impossible 
unless )Cf3, v (v(m)) fits inside the unit box (the macroscopic vessel of the system) 

which is obtained by the scaling of B^v by the factor 1/(2N). Define 

m(P,rj) = min{ m : K,p tTj (v(m)) CD}. 

Theorem 1.1. Assume that (3 £ 2$ and all the components of the boundary magnetic field 
■x] el r have the same sign. Assume, furthermore, that A^ > — Tg(e<j). Then, for every 
m £}m((3,r]),m*((3)[, 

1 Z /3,7? i 
lim — -7—r log — ^' m,+ = lim — ; 1 log uftf 1 , (M»„ < m) 

= -W/3 >v (1Cp tV (v(m))) = -wj (m). 



Theorem 1.1 gives a microscopic justification of the thermodynamical variational prob- 
lem leading to the Winterbottom construction, but, as it is, yields little information on 
the statistical properties of the microscopic spin configuration under the canonical mea- 
sure /i/flj m _|_. In the Li-approach we pursue here the microscopic spin fields are identified 
only on the (renormalized) mesoscopic level through the local order parameters or, equiv- 
alently, through locally averaged magnetization profiles. These local order parameters are 
piece- wise constant functions on the continuum box D and they are constructed from the 
microscopic spin configuration a G { — 1, 1} Djv in the following way: 

Choice of scales. We shall always take the microscopic size of the binary form, N = 2 n . 
Similarly, all the intermediate mesoscopic scales are of the form K = 2 fc , k E N. 

Partition of Bjy. At each fixed mesoscopic scale K = 2 k we split the microscopic ves- 
sel Bjv into the disjoint union of shifts of the mesoscopic box IB a- = {-\K,...,\K} d . 

These shifted boxes are centered at the lattice points from the rescaled set I$n,k = 
K(B N/K -(l/2,...,l/2)): 



where Mk(i) = i + Mr. 



Induced partition of B. We scale (1.12) by the factor N. With each mesoscopic lattice 



box Mx(i) in the partition ( [1.12 ) we associate the continuum box 

K l K~ 



k(x) = x+ 1 
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centered at the point x = jji. We use H>n,k to denote the set of all such centers x; 
^>n,k = jj^n/k- The induced mesoscopic partition of D is given by: 

( - \ 

D = closure \J M NjK (x) ■ (1.13) 

The local magnetization profile Mn,k corresponding to a microscopic spin configuration 
a E {—1, 1}° N is a function on the continuum box B, piece-wise constant with respect to 
the partition ( 1.13j ). For every x £ ^>n,k the value of Mn,k on Mn^k(x) equals to the 



averaged magnetization on the corresponding lattice box Mx(i) with i = Nx; 

M N , K (<r)(y) = a i> for yG 

j€B K (i) 

Informally, Mnk( (T ) is the resolution with which one observes the microscopic spin field 
a. Local proximity of the Gibbs measure to one of the pure phases is, in this way, quan- 
tified by the proximity of Mnk(ct) to the corresponding order parameter ±m*(/3). Our 
probabilistic counterpart of the thermodynamic Theorem |Ll] states that on coarser scales 
local magnetization profiles M n,k comply, with the overwhelming conditional probability 
A 4 !} m + ' w ith the thermodynamic prediction. This agreement has to be understood in 

the Li-sense: Given a measurable subset A C W 1 let us define the function 1^ £ Li(B) 
via: 

+ 1, ifxGA 
-1, iix^A 



1a(x) 



Theorem 1.2. Fix a number v < 1/d and assume that (3 £ 25 and all the components 
of the boundary magnetic field n G R r have the same sign. Assume, furthermore, that 
A^tj ^ —Tpi^d)- Let m £ ]fh(/3, rj), m*{(3)[ . Then for every 5 > fixed, one can choose a 
finite mesoscopic scale Kq = Ko(f3,r],6), such that 

N™oc K ™K% N» ^N,m,+ ( m j n ll-^iV^(^) " ™^ I+ ^,,(»( m ))) \\ U (d) < 5 ) = 1 ■ 

Of course, the above minimum can be restricted to the shifts of K,p„(v(m)) along the wall 
Vd in the partial wetting/drying case; |A/3 i7? | < ra(e*rf), and, respectively, to all admissible 
shifts of the Wulff shape K.p{v) within D in the case of complete drying; Ap )r) = Tg(ed). 

1.7. Structure of the paper and further remarks. To a large extent the proof of 
Theorems [0| and L2 is a book-keeping exercise based on the renormalization formalism 



we have developed for the nearest-neighbor Ising model in |BIV1 |. In that paper we tried 



to decouple deep model oriented facts from the relatively soft Li techniques. The fact 
that the proof swiftly goes through in the general case we consider here is a dividend of 
such an approach. 



The deep model oriented fact is the validity of Pisztora's coarse graining [Pi| and its 



application to the relaxation properties of the FK measures [CePil]; both have been orig- 



inally developed in the nearest neighbor context, but go through literally without changes 
in any model with translation invariant finite-range ferromagnetic pair interactions. In 
Section Q we briefly recall the corresponding construction and show how to adjust it to 
treat the boundary surface tension. In Section || we set up the machinery of random 
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mesoscopic phase labels and relate it to the context of the theory of functions of bounded 
variation. Since the boundary field introduces a singularity of the surface tension at the 
bottom wall, some additional care (as compared to JBIV1 |) is needed at this stage. Fur- 



thermore, the appropriate FK-representation in the case of negative boundary magnetic 
fields happens to be a conditional one. In particular, the compactness estimates have to 
be modified. 

The proof of Theorem |1.2| (and, on the way, Theorem |1.1| ) is relegated to Section ||. In 
view of the results and adjustments of the preceding sections it closely follows the pattern 
laid down in [BIV1]. Accordingly, our exposition in this section will be rather concise with 



repeated references to the appropriate parts of | BIV1 ] . 



Finally, we would like to point out that though the minimizers of the variational problem 
(WP)^ restricted to the sets PCI are, in general, not known for m < mp jV , the Li- 
approach still could be pushed through to yield meaningful results in the sense that one 
can prove that the mesoscopic configurations concentrate around a (in general unknown) 
set of surface energy minimizers. Such an idea of decoupling between the probabilistic 
analysis of the microscopic model and the investigation of the limiting variational problem 
has been put forward in PV1] and, recently, systematically exploited in [ CePi2f| . Indeed, 



the microscopic derivation of phase coexistence is based upon local decoupling estimates 
and on surgical procedures in small regions localized along the interface of the crystals. 
Therefore, from the probabilist point of view the approach is local and does not rely on the 
global shape of the solutions to the corresponding variational problems. On the technical 



level, however, one should adjust the approximation results of Subsection 4.1 to the case 
of unknown minimizers. The approximation procedure we use here hinges on the fact that 
the minimizers are known to be convex. If the surface tension is singular at the boundary 
(that is in the partial wetting/drying cases), a complete geometric analysis without such 
an information would require a serious detour into the geometric measure theory which, 
from our point of view, would not bring any deeper insights on the microscopic phenomena 
involved in the Winterbottom construction. We refer to [Ce| or to [CePi2] for a detailed 



exposition of the appropriate facts from the theory of functions of bounded variation. 

A more challenging problem is to gain a better understanding of the statistical properties 
of phase boundaries, even on mesoscopic or macroscopic scales. For example, the only 
information the Li-approach yields on the entropically repulsed interface is a localization 
on the macroscopic scale, which is, of course, not such a fantastic result, since we absolutely 
do not control the typical microscopic size of the corresponding fluctuations. 

2. FK Representation 

FK representation is an artificial coupling between the ferromagnetic spin models and 
the, so called, random cluster measures. Its power rests with the fact that the latter are 
in the uniqueness regime and possess uniform decoupling and relaxation properties, even 
when the original spin system is in the phase coexistence regionPJ. The art of the FK 
representation is, therefore, to single out such spin configuration properties which admit 
a reformulation in the FK language. 

2.1. Representation of the bulk states. With a given finite range translation invariant 
pair interaction potential J = {Jij = Ju-j\} we associate the graph (Z rf ,£j) where the 
edge set £j consists of all (unoriented) pairs of vertices (i, j) with J\i-j\ > 0. Each such pair 



3 



Strictly speaking, this is known to be true for all sub-critical temperatures, except possibly for a 



countable subset of the latter, see Subsection 3.2 
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will be called a bond of £j. The set fij = {0, is the sample space for the dependent 
percolation measures associated with J. Given uj £ flj and a bond b = (i , j) € £ j, we say 
that 6 is open if uj(b) = 1. Two sites of Z rf are said to be connected if one can be reached 
from another via a chain of open bonds. Thus, each splits Z rf into the disjoint union 

of maximal connected components, which are called the open clusters of f2. Given a finite 
subset B C Z rf we use cb(u;) to denote the number of different open finite clusters of u> 
which have a non-empty intersection with B. 

We define next the FK measures which correspond to the finite volume spin Gibbs states 
on the boxes Bat: First of all these measures put non-trivial weights on the percolation 

configurations uj G £In,j = {0, 1} £jv > j , where the set Sff j comprises those of the bonds 
b £ £j which intersect Bat. The boundary conditions are specified by a frozen percolation 

configuration tt £ 3 = Oj \ CIn,j- Using the shortcut c* n (uj) = cm n (oj V ir) for the joint 

configuration uj V 7r 6 £j , we define the finite volume FK measure $ *T on f2 jy, j with the 
boundary conditions tt as: 

*#» = i ( n (i-ww) 1 ^^^))" (2.i) 

where, for a bond b = (i , j) E £n,j, we define the corresponding percolation probability 
Pp{b) = 1 - exp(-2/3J| i _ i |). 

The measures $h are FKG partially ordered with respect to the lexicographical order 
of the boundary condition it. Thus, the extremal ones correspond to the free (it = 0) 
and wired (tt = 1) boundary conditions and are denoted as <3?^ f and respectively. 
The corresponding infinite volume (N — > oo) limits and $ /3 ' w always exist, and it 
takes a relatively soft argument |Gri| to show that for all but at most a countable set of 
temperatures [5 S M+ there is a unique infinite volume FK-measure: 

(Al) We assume that the inverse temperature (5 satisfies $^ ,f = $^ ,w . 

Assumption (A) is a combination of (Al) and the condition that there is percolation in 
slabs, see Assumption (A2) below. 

The spontaneous magnetization m*{(3) admits the following expression in FK terms: 

m*{fi) = lim ^ W {0^M C N ) = $^' w (0 <-► oo) . (2.2) 

More generally, the finite volume (spin) Gibbs state ^ N + on { — 1, 1} Bjv (see Subsection |1.3|) 
can be recovered [ESJ from the wired FK-state < J ) A ^ W as follows: Sample a percolation 



configuration u> G &n,j from $^ w . Spins at the sites of B^v which belong to the wired 
component of uj VI (that is connected to the boundary B^) are assigned value +1, whereas 
the remaining open clusters of uj V 1 are painted to ±1 with probability 1/2 each. 

2.2. Representation of the surface tension. The surface tension of the spin model 
( 1l.5| ) can be expressed in the FK language in the following way: Let us split the boundary 
dA.M,L{n) = % d \ Am,l(«) into two pieces: 

dA M>L (n) = d + A M:L (n)\Jd-A M:L (n) 

f dA M>L (ri) n {i G Z d : i d > 0}) [j (dk M , L {n) f){ieZ d : 0}) 
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Then the ratio of partition functions in ( |1.5| ) is equal to 

= * 9-A M An)) , (2.3) 

Am,l(") 

where the FK-percolation event {<9 + Am,l(^) *h d~ \M,L(n)} means that no site of Am,l(w) 
is connected both to (9 + Am,l(^) and to 8~Am L(n)- 



The expression ( |2.3j ) for the surface tension is not very convenient in practice. Indeed, 
in higher dimensions one should not expect the interface to decouple along pure wired b.c. 
It would therefore be very useful to have a more robust definition which would still fit in 
the framework of the Li-theory. The idea is twofold. First, one can relax the pinning of 
the "microscopic interface" along the lateral sides of A^l^), by requesting only that the 
top and bottom faces of the box are disconnected. Of course in order to still recover the 
surface tension in the direction orthogonal to re, one has to choose M <C L. The control of 
the interactions on the lateral sides is one of the major technical contribution initiated by 



Cerf in Cel. The second point is to relax boundary conditions, i.e. to replace the wired 
b.c. by some arbitrary tt. In order to avoid a possible deformation of the "microscopic 
interface" due to the attraction by some b.c. tt, one has to let enough room for the system 
to relax to its unique equilibrium phase. This is done by imposing the boundary condition 
tt outside the bigger box A^. 

More precisely, let d top AM,L(n), resp. d hot AM,L(n), be the face of 9 + Am,l(r), resp. 
d~AM,L(n), normal to re. As a consequence of the relaxation properties of FK measures 
derived by Cerf and Pisztora |CePil] (Proposition 3.1), the following holds 



Lemma 2.1. Let e > 0. For any € 25, 

r M = -^J0 1 °Z ^Mft) *4 d bot A eL , L (n)) + c e , L (n, re) , 

with c>e,l (tTj ft) going to zero uniformly in tt and n, as L goes to infinity and e goes to zero. 

2.3. Representation of the boundary states. The FK-notation for the boundary 
Gibbs states closely follows the full space setup introduced in the previous subsection. 
The relevant graph for the inter-spin interactions is (L d , £j) which is just the restriction 
of (Z d ,£j) to the half-space LA In order to incorporate the boundary magnetic field 
rj = (771, rjr) (see (|1.6| )) we augment this graph with a ghost site q connected to all the 
sites in the first r layers of h d . Thus the edge set for the boundary model is given by 



Similarly, given N G N, we use Cfj j to denote the set of bonds of £j which have a 
non-empty intersection with O^. The sample space for finite volume FK states on D^v is 
given by 

= {o, i} £ W 

Assume that r] (that is rj^ ^ for k = l,...,r). Given a configuration tt £ 

rd \ rd o ^ 

{0, 1} J -'7 V the FK measure 3% 011 ^N,J,r) is defined by 



^^(p) = — !— 



11 (i-p^)) 1 -*^)) 6 r^, (2.4) 
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where c^(£) denotes the number of open finite clusters of £ V7r which intersect Bat and do 
not contain the ghost site g, whereas the percolation probabilities pp(b) are defined exactly 
as in (2.1) for the bonds b G C d and equal to 1 — e~ 2 ^ Vk on the ghost bonds b = with 



id = k, k = 1, . . . ,r. We suppress the super-index 7r for the wired state with the boundary 
conditions % = 1. 

The boundary Gibbs state /^jj + can be reconstructed as follows: sample a bond con- 
figuration £ G Hjv,j,?j from <£jy , and paint with 1 all the clusters of £ connected either to 
L d \ Djv or to g, whereas all the remaining clusters of £ are to be painted into ±1 with 
probability 1/2 each. The corresponding joint bond-spin probability measure is denoted 

byP^. ' 

Remark 2.1. The FK measures $^ 7?,7r corresponding to the non-negative magnetic fields 
r] ^ are the basic ones and we shall use them in all our considerations. 

The representation derived in the case of non-positive magnetic fields rj = (— \rji\, . . . , — \ij r \) 
can be described as follows: Define the FK percolation event 5at C Sjv,J?j as 



3tv 



= {t€E N)3 , v | fl ^L d \Biv}, (2.5) 



and set ( • ) = ( - | 3iv) - Then, the boundary Gibbs measure fj,^ 71 + can be 

reconstructed from ( • ) as above, except that the spins connected to q are, this time, 
painted into —1. The corresponding joint measure is denoted by Pjv±- With a slight 
abuse of notation we shall write (keeping in mind that the g-cluster is repainted into —1): 

pftM ( . ) = pAW ( . O . (2 . 6) 



2.4. Representation of the wall free energy. In view of the spin-flip symmetries 
the difference A^n (see ( |1.7|) ) admits the following expression in FK terms (recall the 
convention rj ^ 0): 

A ^ = -^lJ^y^ log ^ v ^ Nh (2 ' 7) 

and, accordingly, if r/ is a negative boundary field A^ = — i^i. 

As in the case of the surface tension, the above definition is too restrictive to be used 
in practice, because of the pure wired b.c. Fortunately, a more robust expression is also 
available here. 

Let d top ' £ Bat = H>N^{i G L d : id = [eN]}, where [x] denotes the integer part of x G M. d . 
Then, using the relaxation properties of FK measures ( |CePil|l , Proposition 3.1), one can 
prove the following 

Lemma 2.2. Let e > 0. For any (3 G 2? and any n ^ 0, 

A/3,r; = -^p T log^( 0t 4a to ^Bv)+c e , J v(vr), 
with c ej v( 7r ) going to zero uniformly in n, as N goes to infinity and e goes to zero. 
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3. Renormalization 

In order to perform the analysis of phase coexistence in a macroscopic setting, a con- 
venient formalism, the Geometric Measure Theory, is introduced in subsection 3A. The 
embedding of the discrete spin system in the continuum setting relies on the introduction 
of renormalized variables, the mesoscopic phase labels (Subsection |3.2|). 



3.1. Functions of bounded variation. We refer the reader to |EG| ] for an introduction 
to the theory of functions of bounded variation and to [BIV1] for a discussion on its 
relevance in the context of phase coexistence. 

Let O be an open smooth neighborhood of B. On the macroscopic scale the system 
is characterized by the boundary condition g G BV \B, {±1}J and by a ±1 phase 

function u G BV(intB, {±1}). The fact that u{x) = 1 for some x G intB means that 
locally at x the system is in equilibrium in the "+" phase. Thus, u should be interpreted 
as a signed indicator of the regions containing the different phases and the boundary of 
the set {u V g = — 1}, where 

| u(x) if x G intB 
u\/g(x) = < ^ (3.1) 

[ g(x) if x G C\B 

as the interface in the presence of ^-boundary condition. 

It is well known [EG] that «Vj £ BV(C, {±1}) whenever the phase function u G 
BV(intB, {±1}). Now, for any v in BV(0,{±1}), there exists a generalized notion of 
the boundary of {v = —1} called reduced boundary and denoted by d*v. If {v = — 1} 
is a regular set, d*v coincides with the usual boundary dv . Given a boundary condition 
g G BV(C \ B, {±1}) and a phase function u G BV(intB, {±1}) we use d*u to denote the 
reduced boundary of u in the presence of the b.c. g: 

d*u = d*{uVg)nB = d*{uVg)\d*g. (3.2) 

In the setting we are working with here, there are two natural types of boundary conditions 
depending on the sign of the magnetic field either g = g + = 1 in the case of positive 
magnetic fields ij ^ 0, or 

f +1, if x d > 
9{X) = 9±{X) = { -1, if x d ^ (3 - 3) 

in the case of negative magnetic fields rj < 0. 

For any non-negative boundary magnetic field r] and any g G BV \ B, {±1}^ define 
the functional 



Jd* g u\V d Jd;unv d 



Set VV^" 7? ( • ) = Wf3 )V (• | g+) and, accordingly, VV^ : 7? ( • ) = W/3 iJ? (• | g±). As we shall see be- 
low this notation exactly corresponds to the color blind formalism of the FK-representation 
and to the construction of the measures P^'^.' in Q2.6| ). Thus, the Winterbottom functional 
( |1.1| ) can be rewritten (recall that for 7] ^ 0, = — A^ i^i) as 

fW+(n) for 

l W /3,|r,|( n ) ~ A /3,h| for 7? SCO 
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In this way the functional Wa^ is unambiguously defined on BV(intB, {±1}) for every 
constant sign magnetic field n. 

For m £]rh(P , rj) , m* the Winterbottom shape JCp tri (v(m)) fits into B and the func- 
tion l£ (v{m)) £ BV(intro, {±1}) is the stable minimizer of W/3 )J? in the following sense: 

For every 5 > there exists e > such that for every function u E BV(intD, {±1}) with 

f^u(x)dx ^ m/m*p, 

Wp, v (u) ^ w* PtJ] (m) + e =^ mill «- 1 X+K ( „ (m)) < 5. (3.5) 

Indeed, let r^ i?? be the support function of the Winterbottom shape itself. Regardless 
of the sign of the monochromatic magnetic field rj define the relaxed functional: 



3+ 



Being the Wulff type functional (or, in an alternative terminology, being the mixed volume 



[SchJ) W/3,n is lower-semicontinuous with respect to the Li(B)-convergence. Furthermore, 
by the refinement jFMj l of the Brunn-Minkowski inequality |Scb( | , the shifts and dilations 
of ICs t n are the only minimizers of W/j « i n the problems with the corresponding volume 
constraints. Hence, Wf3 }V satisfies the stability property ( |3 . 5| ) . However, W^ i7) VVg iT) 
and their values coincide on the shifts and dilations of the Winterbottom shape As 
a result, the stability property ( |3.5| ) is immediately inherited by the original functional 

In the partially wetting case < T p{^d) the minimum in ( |3.5| ) is over all admissible 
shifts x £ Vd along the substrate wall. In order to see this notice that in the latter 
situation the minimum (m) is strictly less than the unconstrained minimum of the 

Wulff functional Wp{u) over the functions u £ BV(intB, {±1}); J^u(x)dx ^ m/m*Q. 

Fix now a point x £JS>\Vd and define r = Xd > 0. For any 5 > and any function u £ 
BV(intB, {±1}) satisfying \\u — I^-f/c^ v (y(m)) ll^rfi) ^ ^ one choose a section ^ r such 
that 7i 1 ({2;^ = rs} Pi {x : u{x) = —1}) ^ 5/r. Modifying u as us = +1 on < r$ we 
obtain that Wp(us) ^ (tx) + 5 max rg/r and, of course, that f^us(x) dx ^ (m+5)/m*^. 
As a result, a W/j^-minimizing sequence cannot converge to l x+ ic /3r] . 

3.2. Mesoscopic phase labels. In the Li-approach the local order parameter Mn,k is 
quantified by three values { — 1, 0, 1} according to the local proximity of the system to one 
of the pure phase. Such a renormalization procedure is delicate to implement directly. 
One way, proposed by Pisztora in |F]] , relies on the FK representation of Ising model. A 
technical advantage to work with the FK representation, is the uniqueness of the measure 
in the thermodynamic limit even when there is breaking of symmetry for the Ising model. 
In a sense the FK measure is much less sensitive to boundary effects than the Gibbs mea- 
sure, which underlines the required decoupling properties. 



Typical configurations for the FK representation : 

As explained in (|2.2j), the FK counterpart of spontaneous magnetization is a uniform 
positive probability that the site is connected to the boundary of arbitrarily large boxes 
with wired boundary conditions. The implementation of Pisztora's coarse graining requires 
an a priori stronger notion, namely percolation in slabs : 
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(A2) 35 > 0, 3L > 0, VL > L Q , lim^oo wf XtVe s LtN $s'l N ( x v) > 6 > 

where Sl,n is the slab {|ii| ^ L, — N ^ ik ^ N, k = 2, . . . , d}. The importance of this 
notion was first realized in [ACCFR| in the context of Bernoulli percolation. If (A2) is 
satisfied then there is percolation and (3 > (3 C . The critical value (3 C above which (A2) 
holds is called slab percolation threshold. We define 03 as the subset of ]/3 c ,oo[ for which 
assumption (Al) holds. The set 03 differs from ]/3 c , oo[ by at most a countable number of 
points. In fact, it is conjectured that 03 =]/3 c , oof. 

Originally (see W%) the estimates on typical configurations were devised in a nearest 
neighbor set-up. Inspection of the proof shows that the results extend readily to models 
with finite range interaction. 

Property (A2) enables to obtain enhanced estimates on the percolation in domains 
by slicing these domains into slabs. Furthermore assumption (Al) implies a relaxation 
property of the system and therefore uniform estimates can be deduced for arbitrary 
boundary conditions. 

We can now proceed to describe Pisztora's coarse graining. Let a be in (0, 1) and fix 
Q > 0. On a coarse grained scale K € N, the typical configurations u in the box Mk satisfy 
the 3 properties below : 

1. There is a unique crossing cluster C* in Mk, i-e. a cluster which is connected to all 
the faces of the vertex boundary 

dM K = {x G M C K | 3y G Bjc, such that J(x - y) > 0}. 

The faces are simply the faces of the box Mk enlarged by a factor depending on the 
range of the interaction. 

2. Every open path in Mk with diameter larger than K a is contained in C*. (The 
diameter of a subset A of Z d is sup x yeA \\x — y||i-) 

3. The density of the crossing cluster in Mk/2 is close to m* with accuracy £ > 

|c*|€K-c,m* + c]fr, 

where | • | denotes the number of vertices in a set. 
A configuration u) in Mk is good if it satisfies the 3 assertions above. It was proven in 



[Pi| that for any € B and for large enough mesoscopic scales, good configurations are, 
uniformly over boundary conditions, typical in the following sense : 

inf <&j? (u is a good configuration in B^) ^ 1 — exp(— crK a ) , (3-6) 

7T 

where > depends on the accuracy £ of the coarse graining. 

Mesoscopic phase labels : 

On the mesoscopic scale K = 2 k , each box Bat i ^(x), centered at x & ^n,k is labeled by 
the variable Ujf^fajX) 

cfk, \A |l if oj is a good configuration in ~B>2k(Nx), 
G Djv,k, uy fK (u),x) = < 

I otherwise, 

where the information recorded in u^ K (u,x) takes into account the bond configuration 
in the larger box M2k(Nx). We say that a block Mn : k(x) is regular if u^^x) = 1. 
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Estimate ( |3.6| ) is uniform over the boundary conditions, therefore it enables to control 
the variable u^f K (x) independently of the events which occurs on boxes which are not 
*-connected to Bjv,jf(sc). Following [jpj, RLSS |, this implies stochastic domination of the 



the field {x G H>n,k I ^^(a;, sc) = 0} by the Bernoulli site percolation measure with 
parameter exp(— c^K a ) on On,k', in particular, 



*N W '" ( u n F k(^) = °' • • • >«&f* l*t) = 0) < exp(-c c K a £) . (3.7) 



Furthermore, the basic FK-estimate ( |3.6| ) enables to control the local magnetization 
profiles Mn,k, f° r large enough mesoscopic scale K = 2 k , by a coarse graining which 
satisfies similar decoupling properties. Let x be in B>n,k- Then for any y G Mn,k(x), we 
define the mesoscopic phase label u N K as: 

C ( n A J sign(C*) if Uj^ K {u,x) = 1 and I-M^kCo", x) - sign(C*)m*| < 2Q, 
u NiK {a,u,y) - ^ otherwis e, 

(3.8) 

where C* is the crossing cluster in B^(x). 

A block Mn,k(%) may have label for two reasons : either K x) = 0, or a; is a good 
configuration on ~M>2k{Nx), but the local average magnetization M.n,k{%) on B^(A r x) is 
non typical. The former case is covered by ( |3.7| ). In the latter case, the corresponding 
shift of the magnetization is entirely due to an abnormal coloring of the small (i.e. with 
the diameter less than K a ) w-clusters intersecting Mk(Nx). The probability of such a 
deviation is, for each particular block Mk(Nx), bounded above by exp{— c( 2 K d ( 1 ~ a ^}. On 
the other hand, all small clusters intersecting Mk(Nx) lie inside M2k(Nx), which already 
leads to the required decoupling properties. Notice that the coloring argument inside a 
good block is insensitive to the repainting of the g-cluster into —1 on the event 3/v- We 
refer to the original article [|P| for the precise workout of these estimates and to [Bl] where 
the proof has been adapted to the setup we employ here. 

As a result we obtain that there there exists a sequence {p^} with ]imK-> o0 p* K = 0, 

such that the field {x G H>n,k \ u^ N K (a, u, x) = 0} is stochastically dominated by the 

'c 

Bernoulli site percolation measure Pp^i-c on B>n,k- I n particular, for rj ^ 



■ N 



(u^ K ( Xl ) = 0,...,u^ K (x e ) = 0) ^(piY. (3.9) 



Remark 3.1. If for two different points x,y G D>n,k the corresponding u N K -phase labels 
have opposite signs, that is if u^ N K (x)v!j^ xiv) = — 1> then on any finer scale K' ^ K any 
^-connected chain oflS>N,K' blocks joining Mn^k(x) to Mn,k(u) contains at least one block 
with zero K , ( and, hence, with zero u^ N K ,)-label. This follows from the fact that 2 boxes 
with opposite labels cannot have a common crossing FK- cluster; thus they are necessarily 
separated by a contour of 0-blocks. 

Remark 3.2. The modification of the boundary interaction has no impact on the coarse 
graining. Since the slab percolation threshold hypothesis requires percolation on slab with 
free boundary condition, the presence of a positive boundary magnetic field can only im- 
prove the estimates in view of FKG monotonicity. Furthermore, in each box touching 
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the wall, the event that a crossing cluster is connected to the ghost point q occurs with 
overwhelming probability (the slab percolation hypothesis (A2) is valid uniformly over the 
boundary conditions). Finally the control of the density of the crossing cluster (property 3 
in the definition of a good block) is unchanged because the boundary field has no influence 
on the bulk properties due to the hypothesis (Al). 

The mesoscopic phase labels u* N K are {±1, 0}-valued function on D which are piece- 
wise constant with respect to the partition ( [1.12 ). We are now going to prove that the 



Li-difference between the local magnetization profiles Mn,k and the phase labels u* N K 
vanishes on the exponential scale. Thus, from the point of view of the Li-theory these 
objects are indistinguishable, and one is entitled to switch freely from one to the other. 
In particular, the fundamental exponential tightness result which we shall establish below 
in Subsection |3.4] comes naturally in the context of the {±1, 0}-valued phase labels. 

3.3. Relation to magnetization profiles. 

Lemma 3.1. Let a positive boundary magnetic field rj ^ and a number v < 1/d be fixed. 
For any 5 > 0, one can choose the accuracy C of the coarse graining and a finite scale 
Kq = Ko(5,f3), such that 

iv^/.oW^ 1 ^^^ 11 ^^"" 1 ^ 111 ^^ 5 ) ^ -<CK)N l - d ». (3.10) 

Remark 3.3. By ( p. 7] ) and ( |2.6| ) the super surface order exponential estimate ( |3.10| ) 
holds, in the case of negative boundary magnetic field rj ^ 0, also for the conditional 
measure IP^' ± . 

Proof. Recall that the coarse graining is initially defined in terms of the mesoscopic 
FK-variables u^^(uj) which depend only on the bond configurations. The phase label 

u^ N K (a, uj) takes also into account the random coloring of good FK-blocks. In fact, 



\\M N ,K-m*u%, K h ^C + - Yl 1 



zero W jc 



and the claim of the lemma follows by the domination by Bernoulli site percolation (3.9), 



actually for any choice of £ < 5. □ 

3.4. Tightness theorem for mesoscopic phase labels. We shall formulate a tightness 
result which holds under fairly general phase boundary conditions outside D. The proof 
is a straightforward adaptation of the argument used to verify the claim of the tightness 



Theorem 2.2.1 in [BIV1] and, as we shall briefly indicate in Remark below, it can 



be trivially extended to cover the case of g-valued phase labels as, for example, recently 
considered in [CePi2|. 

Given N = 2 n and K = 2 k (n ^ k) consider, as in ( |1.12j ), the following splitting of M. d : 

V ®N,k(x)> (3-11) 



In what follows we shall use the same notation as in Subsection 3.1. Thus, O is a fixed 
smooth neighborhood of D. A boundary condition is an element g G BV (o \ D, {±1}^ , 
which is piece-wise constant with respect to the partition ( 3.1 1] ) on some (M, L)-scale. 
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Notice that such boundary condition g is automatically piece-wise constant on all further 
(N, K) scales with N^M and N/K > M/L. 

A family {un,k} of mesoscopic phase labels is said to be compatible with a boundary 
condition g; un,k ~ 9, if 9 is piece-wise constant on the (iV, if)-scale and the following 
phase rigidity condition (R) is satisfied: 

(R) Define wn^k = un^k V g as in fl3,l| ). Then for every x € ~Bn,k and y G 2^ if 

with wn,k{x)wn,k{u) = —1 and for every L ^ K any *-connected chain of Bjv.L-blocks 
leading from Bjv,A"(ic) to Mn,k(u) necessarily contains at least one B^L-block with the 
UAr.i-label zero. 

In order to state the tightness result we need to define the perimeter of a BV-function 
u £ BV(intB, {±1}) under a boundary condition g: 

Vg{u) = U d -l (<9>) • 



It is well known [EG| that for every a ^ and for every boundary condition g the set 
C g (a) = |u G BV(intD,{±l}) V g (u) < a} 

is compact in Li(D). Let P 9 be the minimal perimeter of the functions compatible with 
the boundary condition g. 

Before deriving an exponential tightness for the mesoscopic phase labels under measures 
with boundary fields, we first state a general result for an abstract measure P and general 
boundary conditions 

Theorem 3.1. Let a sequence of Tiofi-Ticgativc fiuTnbcrs satisfy \iiB.^^ 00 

Pk = 0. 

Assume that a family of random mesoscopic phase labels {itjv k} on ^ satisfies the follow- 
ing decoupling condition (D): 

(D) For every N and K the zero-label field |x £ K | uj\r t K( x ) = 0j under the mea- 
sure P is stochastically dominated by the Bernoulli site percolation measure Pp^c- 

Fix a number v < 1 jd. Then for every 5 > and for each boundary condition g there 
exists a finite scale Kq = Ko(d,5) and a constant c = c{d) > 0, such that, uniformly in 
a > P g , 

1 c 
limsup max — ^ logP [u N>K ~ 9 \ u N)K V (C 9 (a), 5)] < - — j^a , 

where V(C g (a),S) denotes the 5 -neighborhood in ofC g (a). 

Remark 3.4. TTie tightness result above could be trivially extended to the case o/{0, 1, g}- 
valued phase labels. One needs only to modify the rigidity condition (R) (to require a zero 
block in any ^-connected chain between two blocks with different non-zero labels), and rede- 
fineVgiu) as the TL^ d ~ 1 ^ -measure of the jump set ofu\/g [AmBr]. Ifujsr t K is a {0, 1..., q}- 



valued phase label and g is a boundary condition which satisfies such modified assumption 
(R) and the decoupling assumption (D), then, for every i = l,...,q, the {0, ±l}-label u l N K 
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defined as 

if un,k{x) = i 
u n,k( x ) = { if un,k(x) = 

otherwise 

and the boundary condition g % similarly defined already satisfies the assumptions of Theo- 
rem \3.J{ . It remains to notice that 

i 

{u N , K e V (C g (a),8)} D p| {u% )K e V (C gi (a/q) , 8/q) } . 

i=l 



Proof. The proof essentially follows Theorem 2.2.1 in [BIV1]. The only difference with the 



periodic case discussed there is the existence of open contours attached to the boundary. 
But, for any choice of the macroscopic boundary condition g (as defined above), these 
boundary contours are, in the language of [BIV1], large, that is their diameter exceeds 



Cd log N on all sufficiently big (N, K) renormalization scales. Let us briefly recall the 
three main steps of the argument: 

The first step is to check that the volume of 0-blocks is negligible. This simply follows 
by the assumption (D) : 

p|#{i6D v :b^(i) = 0}^^J J < ex P^(^) log (^)) • (3 ' 12) 

The second and the third steps are devoted to the control of the regions surrounded by 
contours, i.e. by connected surfaces of 0-blocks. By phase rigidity assumption (R) the 
contours play the role of mesoscopic interfaces separating regions with mesoscopic phase 
labels of different signs so that they contribute to the total perimeter of the configurations 
Un,k- We distinguish between two types of contours, the small contours with diameter 
smaller than Cd log N (where Cd is a given constant) and the remaining contours, namely 
the large ones. 

Peierls type estimates can be applied to bound the probability of events such that the 
collection {Tj} of large contours has a total area larger than a{^) d ~ l 

E i r *i >(§) a ) < ex p \-^ a {jc ) ) ■ ( 3 - 13 ) 

. Y{ large / \ / 

The small contours are not controlled in terms of their total area but of the volume of 
the regions they surround. By the choice of the macroscopic boundary condition g there 
are no open small contours, and, proceeding exactly as in the proof of Theorem 2.2.1 in 



[BIV1], one can show that 

E voLCTO^) ^exp(-c 2 ^^). (3.14) 

\T % small / 



As a consequence of the previous estimate, the small contours will have no contribution 
to the Li-norm. Combining estimates ( 3. 12|) , ( |3.13j ) and ( 3.14| ), we arrive to the claim of 



the Theorem. □ 
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Remark 3.5. The last step of the proof should be modified for general, that is not nec- 
essarily piece-wise constant with respect to the partition Q3.ll ), boundary conditions g G 



BV [O \ B, {±1}J . See also [|CePil[| and [|CePi2f| where small contours have been ignored 
and, accordingly tightness estimates have been derived only on diverging (N, K)-scales. 

Let us go back to the phase labels {m 1 ^ k } we consider here. The boundary condition 
g is given by g = g + = 1 in the case of positive boundary magnetic fields rj ^ and by 
g = g specified in (^) in the case of negative rj < 0. 

Proposition 3.1. Let (3 E 58. Then there exists C{(3,d) > such that for all 5 positive, 
one can find Kq = Kq(5, f3,d), £ = C(5) such that 

1. In the case of positive magnetic fields i] ^ 0, 



lim sup 



1 



max 



log 



■ N 



u% K ^V{C g+ {a),5) 



2. In the case of negative magnetic fields n < 0, 



C 



K, 



d~l 



a. 



limsu i m » 



logP 



AN 

JV,± 



K ?V(C g ±{a),6) 



d-l 



a + |A 



Proo/. Let us briefly work out the rj < case: Recall that pJ)J( ■ ) = pj ,7?l ( • ; J^) /P^ 1 ( 3 at 
On the other hand, any phase label ifi N K which is compatible with 3^ is, after the repaint- 
ing of the g-cluster into —1, automatically compatible with the boundary condition g 
and, therefore, satisfies the conditions of Theorem [O] (under the measure ffjy '). Thus, 
the claim of the proposition follows from the general tightness Theorem [O] and, in the 
case of negative magnetic fields n < 0, by the representation formula for the boundary 
surface tension fl2.7|) , □ 



4. The proofs 

We going are to derive sharp asymptotics involving the surface tension. As mentioned 



in Subsection 3.4, the case of negative boundary field requires more care, therefore we 
focus on this case in the following. 

4.1. Approximation result I. Given a boundary magnetic field rj and m e]mg jr) , m^[, 
we are going to approximate the Winterbottom shape K,p^(v(m)) by regular sets (see 
figure ||). 

Lemma 4.1. For any 5 > one can construct a polyhedral set fCp satisfying: 



and 



< s. 



Proof. It is a standard result of the theory of convex bodies (see e.g. |Sch| ] Theorem 1.8.13) 
that for every 5 > there exists a convex polyhedral set Kg such that the Hausdorff 
distance dnausd^rt; fcp) ^ 8. Of course, the support function ri of such /Ci satisfies 



max 



s$ 5. 



(4.1) 



Define now 



ICi n {x : x d ^ - A^}. 
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Figure 2. Left: The approximate Winterbottom shape K,p „(y(m)) and the boxes of 
Proposition 4.1. Right: A function v G BV(intD, {±1}) (the shaded area corresponds to 



1), and the boxes of Theorem |0] . 

A comparison with reveals that 

dHa US d(^, ?? ,^/3,r ? ) O and max |r| Jn) - Tp >r) (n)\ < S, 



(4.2) 



where and Tp )T) are the support functions of Kffe and JCp >r} respectively. Now, for 
every bounded convex K, C M. d the volume vol^ (JC) can be written in terms of its support 
function tjc as 

vobiC/C) = i / r^n^dWf 1 '. 



Consequently, (4.2) implies that 



€ d5\d)C 



P,v\ 



As we have already mentioned in the end of Subsection |3.1| , 

T PtV (fi x )d'H<t 1) = Wp, v (Kfi, v )- 



die 



On the other hand, for n = — e<2, the support function r| equals to A^. Moreover, for 
every x G <9/C^ \ {y : yd = — A^} any supporting hyperplane to dJC^ at x is also a 
supporting hyperplane to cWCjg. It follows that the support function (n x ) = T^(n x ) on 



x E dfcpn \ iv '■ Ud = — A^}. Thus, as it follows now from (|P| 



< 5 



dIC 



3-V 



Since by (4.2), |vold(/C^ ) — voL^/C^)] ^ <5 1 cJ/C^g^ j , and both |9/C/3^| and \dJCi | are 
bounded above by some finite constant ci(/3, d) < oo, we finally obtain the following 
estimate 

Wp, v {K s p ) - Wp,r,{1Cp, v ) < c 2 ((3, d)S, 
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and the claim of the theorem follows once we redefine /C^ = /C^ C2V1) . □ 
4.2. The lower bound. As usual, we fix a number < v < 1/d. 

Proposition 4.1. Assume that (3 £ 2$ and all the components of the boundary magnetic 
field 7] £ M r have the same sign. Assume, furthermore, that > — Tp(ed). For every 
5 > 0, there is a finite scale Kq{5) such that for every m e]m(/3, rf), m*{j3)[, 

1 

^N^oo N d ~ l K (S) <C N" 

where minimal surface energy value w^^m) has been defined in ( |1.11| ), and the function 
o(-) depends only on (5 and rj and vanishes as the resolution 5 goes to zero. 

Proof. We shall give the proof only in the more difficult case of negative boundary fields 



at5=i v w log ^I+ W Ml f>K ~ ^^KtnJjIlLi < 6) ^ - w}, n (m) - o{5) 



rj < 0. Also according to Lemma 3J-, it will be sufficient to derive the proposition with 
■Mn,K replaced by m*u^ N K (for ( small enough). 

Starting from the approximate shapes Kp let us use the transformation ( |1.3| ) to define 
the scaled polyhedral approximation JCp (v(m)) of JCp iV (v(m)). By the approximation 
result of Lemma [4.1| it will be enough to prove Proposition |4.l| with (v(m)) instead 
of JCp )Jt (v(m)). 

Define = Va nD. Since JCi (v(m)) is polyhedral and convex it is an easy matter to 
show that there exists a side length h > and a finite number t of disjoint parallelepipeds 
R},...,R e with bases B 1 ,...,^ included in d ± tC 5 fjri = dK S p v A 7?f of side length h and 
height Sh such that: 

la) The sets B 1 ,...,^ cover d^Ki up to a set of measure less than 6 denoted by 
U 5 = d ± K,^ v \ U - = i & and the y sat i sf y 

| rtfiiWl + A^i 1^*1 ~ W AN <5 > (4 - 3) 

where n, is the unit normal to B % (and, hence, to the corresponding facet of d^lCp ). 
Ib) For any i = 1,...,£ the base B l divides R 1 into two parallelepipeds R l ' + and 
such that in the case of Hi / -e d , R l >- C K s p v and R^ + C 0\/C^, whereas in the case of 

Hi = —ed the corresponding base B % C V^, while R h+ C D\/C]| and C {x : < 0}. 
In order to enforce a microscopic interface close to the polyhedral set Kp , we define 

i 

A=f]{d-R l N ^d + R t N }. (4.4) 
i=l 

Let us also introduce the set A chosen such that for any configuration £ in A Pi *4 the 

'V 

^={^6 3^, e& = 0, if dist(&,L#)<10}. 

We get 



bonds inside and outside NICp are decoupled. For this, it is enough to close the bonds 
which are in a neig hborhood of U S N = NU 6 (the microscopic counterpart of U 



^t(h C N,K ~ ^.Jk < 6 ) > < l ,l(W<K - I^Jk < ^^i) -(4.5) 
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We first check that 



5 

liminf min P^± ~ Itf Ik < o I ^ n > 3 ■ ( 46 ) 



From the argument developed in the proof of Proposition |3.1| , we are going to check 
that in each of the two connected components of D \ /Cg the phase labels u^ N K are, with 
a uniformly positive probability, close to the equilibrium values (for the Li-norm). 

Since AC\A decouples the connected components, it is enough to consider the mesoscopic 
phase labels in the interior of the regions decoupled by A D A. As the boundary of 
d^lCp is regular, these regions are fairly exhausted by the mesoscopic boxes for any 
scale K 6 [Kq,N v \. Using the terminology of Proposition |3.1| , we assert that with large 
probability, the mesoscopic label configurations are almost uniformly constant and contain 
only contours which have no contribution in terms of the Li-norm. According to estimate 



( 3.14 ), the volume of the regions surrounded by the small contours is arbitrarily small in 



the thermodynamic limit. Furthermore, from the usual Peierls estimates, the volume of 



the large contours is negligible as well. This implies (4.6) 
Finally as A D A C 3, we see by applying (|4.6|) that 



The denominator can be easily estimated by ( |2.7l ) 

d>J |r?l ' w p) = exp ( - A m N^ 1 + oiN*- 1 )) . (4.8) 

On the other hand, as we are trying to bound probabilities from below, even under the 
conditioning on A one can still apply FKG arguments (non-crossing is a decreasing event) 
and use fl2,3| ) and ( p.7[ ) to decouple between different I?j-events which constitute A in 
(|4.4| ). Therefore, since, $^' ?? '' w (*4) ^ e -c3SN d 1 ^ we are en titled to conclude that 



$^ r,lw {AnA)> expl-N^i rtfW&l + A m £ \B l \) - (c 3 6 + o(6))N d - 1 

(4.9) 



Combining fl4.9| ) with the estimates Q4.7| ) and Q4.8| ) ; the approximation property (fD^) ; the 
choice of the polyhedral set ]Cp (v(m)) through Lemma 4.1 and the representation (|37 



of we arrive to the claim of Proposition |4.1. □ 



4.3. Approximation result II. Approximation results used in the proof of the upper 
bound should apply to any v £ BV(intB, {±1}). On the other hand, on the level of the 
precision provided by the Li-theory the underlying probabilistic estimates are less delicate 
than those needed for the proof of the lower bound for the distinguished almost optimal 
polyhedral shapes of the preceding subsection. In fact, global bulk relaxation properties 
play no role at all, and the Li-upper bound is just a coarse estimate on the surface 
tension price of the localized interface. In particular, the approximation construction of 
|B1] suffices and the proof of the corresponding result in the latter paper literally goes 



through. As in the case of the lower bound let us concentrate on the more difficult case 
of the negative boundary magnetic field rj 0. We claim: 
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Theorem 4.1. Let v G BV(intB, {±1}). For any 5 positive, there exists h positive and a 
finite number £ of disjoint parallelepipeds R 1 , . . . , R with basis B l , . . . , B^ of size h and 
height 5h such that: 

II a) For every i = 1,...,£ the base B l divides R l into two parallelepipeds R l ' + and R l ~ 
and we denote by Hi the normal to B l . Either B l n <9B = and R 1 is included in B ; or 
B l is included in 5B \ Vd and only R % '~ is included in B, or B l is included in Vd and only 
R l ' + is included in B 

lib) The parallelepipeds R 1 ,... ,R approximate the reduced boundary d*±v (see Subsec- 
tion 



3.1) in the following sense: 



X ni (x) - v V g ± {x)\ dx ^ 5vol(ff ), 



R 

R l 



where X^ = 1^,+ — and the volume of R l is vol(i?*) = 8h . Furthermore, 



As we have mentioned, Theorem |4.1| (see figure ||) is proved as in the standard (+ b.c.) 
case, see pl[ |. 

4.4. The upper bound. 

Proposition 4.2. We fix (3 G 23 and r] < 0. For all v in BV(intB, {±1}) such that 
Wp,r)(v) is finite, one can choose 5q = 5q(v), such that uniformly in 5 < 5q 

lim sup — - max \og^\ {\\M N ,K - m*v\\ u ^ 6) ^ - WpJv) + o(5) , 

N^oo N d ~ l K (S) ^.N" ,+ 

where the function o(-) depends only on (3 and v and vanishes as 5 goes to 0. 

We fix if large enough and (5 > 0. Using the approximation Theorem, there is 5o < S 
such that for any 6' < So 



< limsup^log/^ G V(#,2*vul# 



= limsup-^ T logPg ^ G nV(^,25vol^) 
where V(R\e) is the e-neighborhood of X^ 



R< 

V{R\e) = tv' G Li(B) | J \v'(x) - X^{x)\ <lv < : \ . 
According to Lemma |3.1|, there exist Kq and £ such that 



lim sup — j-j max logPW — G P| y(if 2<5voli2 i ) 

< limsup JLj Ko max ^ logP^ ^ G f| V(# , 35volff ) 



i=l 
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The first step to extract the surface tension factor is to localize in each R l N a mesoscopic 
interface (where R l N is the microscopic counterpart of R % ). This is done by means of a 
surgical procedure : in each R l N , the Li-constraint u^ N K S V(R l , 3<5vol(i? 2 )) enables to find 

two sections in R 1 ^ and in R 1 ^ (the minimal sections) containing only a small portion of 
the mesoscopic interface. So that up to a small cost one can rearrange the configurations 
in these sections in order to identify clearly the location of the interface. At that point 



the microscopic structure of Pisztora's coarse graining becomes effective (see Remark 3.1 ) 
and the presence of a mesoscopic interface is enough to ensure that the sets d top R l N and 
are disconnected on the microscopic level, where R 1 ' is the parallelepiped included 
in R l with basis B l and height Therefore, introducing the set 

A = jw S 3 : 3a such that u C N)K (a,u) G Q V(R\ 2,5vo\{R 1 ))^ , 



we obtain 



*Sf2 [f]{d top R%^d bot RU) , (4.10) 



The constant C in the error term is proportional to the perimeter of v. It is important to 
note that the surgical procedure is insensitive to the boundary effects and applies equally 
well near the wall Vd or in the bulk. 

Equation (4.10) can then be bounded by 

<&5 W ' W (a{d top i4 t4 <9 bot <}) 



(4.11) 



As in the proof of the lower bound, the denominator can be estimated by (2.7). In order 
to recover the surface tension in each box R l N , we fix the boundary conditions outside each 
box R l N . At that point the relaxed definition of surface tension is crucial and we apply 
Lemma |2.1| or 2.2 depending whether the box R l N lies on the wall or in the bulk. This 



provides the following upper bound 



AM, 

N 



(A) < exp < 



-N 



d-l 



T P (ni)\&\ +A m \ & 



C5N a 



(4.12) 



Using the approximation Theorem ^lj, we, thereby, obtain: 

$^> V (A) < exp{ - N^Wf^v) + o(6)N d - 1 } . 
In view of ( |4.8| ) and the expression ( |3.4|) for the functional VV^i^i the proof of Proposi- 



(4.13) 



tion 4.2 is completed. 



□ 



4.5. Proof of Theorems 1.1 and 



The Theorems 1.1 and |1.2| follow from Proposi- 



tions |3.1| , 4.1, 4.2. We will focus on the derivation of Theorem |1.2| in the case of negative 
boundary field, the other cases can be deduced easily. 
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The first step is to derive the following lower bound for m G]m(/3, rj), m 



Ijmjnf log /zj'^ + (M Djv < m) > - w^(m) , (4.14) 

where (m) was introduced in ( 1.11| ). 

Fix e in ]0, mp — fh((3, rf)[. For any 5' > small enough and a in (0, one has 



On, 



(M 0n < m) > (||A^Ar,AT« - ^I^Wm-eWllLi < 



By using Proposition 4.1 and letting 5' go to 0, we see that 

iimmf -JL- log (M DiV < m) > - w^(m - e) . 

As the minimizers K,p„(v(m)) are obtained by dilatation of the same unnormalized Win- 
terbottom shape tCp^, it follows that 

J™ w *p, n ( m ~ £ ) = w p, v ( m ) ■ 

This concludes the proof of ( 4.14| ). 

In order to derive the upper bound, we define for a given 5 > 

F a = jueL^B) inf ||u — Ix+jc A ,(«(m)) lli-i > *J f~ u ( y "> dy ^ ' 
For a given a > 0, 5' > 0, 

Al + {^ M N,K € F,) *S /igj >+ (^-M^ G F, n V (C 9± (a),5')) 

+/43, + (^^^V(C ff± (a),*0) • 



According to Proposition 3.1, there is Kq depending on 6' such that uniformly in a 
limsup 1 max logP^f^^0V(C 9± (a),y)) < - J^_ a + IAft.,1 



We fix e > and a much larger than £ (wp^m) + 1 A^jJ). The exponential tightness 
will enable us to consider only a finite number of macroscopic configurations: 

Since C g± (a) n F5 is compact in Li(D), we, using the upper bound Proposition 4.2 , can 
cover it (and hence also V (C g ±{a) nF$,5'^ for 5' sufficiently small) with a finite net of 
neighborhoods (V(ui,£i)). <e and choose a finite scale K\{e) such that for every i ^ I 

limsup — I max log fi^f, ( — Mn,K € V(«i, £*) J < - Wa Jitj) + £. 

By the stability property ( |3.5| ) . 

inf W^Cit) > (m) . 

Therefore, by choosing e small enough, we can find a mesoscopic scale K2 and a positive 
constant c = c(5) > 0, such that 

The above inequality combined with the lower bound ( |4.14j ) leads to the result. □ 
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